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ABSTRACT 

Galaxies with the mass of the Milky Way dominate the stellar mass density of the Universe but it is uncertain 
how and when they were assembled. Here we study progenitors of these galaxies out to z ~ 2.5, using data 
from the 3D-HST and CANCELS Treasury surveys. We find that galaxies with present-day stellar masses 
of log(M) « 10.7 built ~ 90 % of their stellar mass since z ^ 2.5, with most of the star formation occurring 
before z ~ 1. In marked contrast to the assembly history of massive elliptical galaxies, the centers and outer 
parts of the galaxies built up at roughly the same rate between z ^^ 2.5 and z '-^ 1. We therefore conclude that 
a "standard" model for the formation of spiral galaxies, with the bulge assembling first and the disk building 
around it, is probably not correct. Instead, bulges (and black holes) likely formed in lockstep with disks, 
through bar instabilities, clump migration, or other processes. We find that after z '^ 1 the growth in the central 
regions gradually stopped and the disk continued to build, consistent with recent studies of the gas distributions 
in z '^ 1 galaxies and the properties of many spiral galaxies today. 
Subject headings: cosmology: observations — galaxies: evolution — Galaxy: structure — Galaxy: formation 



1. INTRODUCTION 

The Milky Way is a very typical galaxy, in the sense that a 
randomly chosen star in the Universe is most often found in 
a bulge - disk system of similar mass. Despite their ubiquity, 
and our exquisite knowledge of one example of their class, 
the assembly history of large spiral galaxies is still very un- 
certain (see Rix & Bovy 2013, and references therein). Disks 
are thought to form through gas accretion, with the structure 
of the disks largely determined by the angular momentum of 
the dark matter halo (Mo, Mao, & White 1998). In many 
models for galaxy formation bulges arise from the merger of 
two disks, after which a new disk can form (e.g., Kauffmann, 
White, & Guiderdoni 1993). In the context of these models 
the Milky Way bulge formed before the disk (see also, e.g., 
Zoccali et al. 2006; Okamoto 2013). This "inside-out" sce- 
nario is qualitatively similar to the assembly history of mas- 
sive ellipticals, which formed a dense core at high redshift and 
then built up their outer parts (e.g., van Dokkum et al. 2010; 
Hilz, Naab, & Ostriker 2013). 

However, recent theoretical and observational studies have 
challenged this picture. In cosmological simulations of gas 
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accretion the structure of the forming galaxy not only de- 
pends on the properties of the dark matter but also on the 
details of the feedback mechanism (e.g., Agertz, Teyssier, 
& Moore 2011; Brooks et al. 2011) and on the accretion 
mode (e.g.. Sales et al. 2012). Furthermore, major mergers 
are probably too rare to form many bulges (e.g., Kitzbichler 
& White 2008), and several studies have suggested alterna- 
tive ways to build up central mass concentrations. In partic- 
ular, bulges may be built by "direct injection" of gas in cold 
streams (e.g.. Sales et al. 2012) and/or by the migration of 
star forming clumps in unstable disks (Elmegreen, Bournaud, 
& Elmegreen 2008; Dekel et al. 2009; Krumholz & Dekel 
2010). Such clumpy, unstable, rapidly star-forming disks have 
been shown to exist at high redshift (e.g., Genzel et al. 2008; 
Forster Schreiber et al. 201 1). 

In this paper (and in Patel et al. 2013b, which uses a com- 
plementary approach) we provide new constraints on the as- 
sembly of disks and bulges by studying plausible progenitors 
of Milky Way-mass galaxies in the 3D-HST survey (Bram- 
mer et al. 2012). The immediate goals are to determine the 
typical star formation histories of these galaxies and to assess 
when their major structural components were assembled. The 
data also provide key constraints on the ingredients in recent 
hydrodynamical models: these models now succeed in repro- 
ducing many of the properties of the present-day Milky Way 
(Brooks et al. 201 1; Guedes et al. 201 1) and to improve them 
further we need to test their predictions at earlier times. A 
Kroupa (2001) IMF is assumed throughout the paper. 

2. MASS EVOLUTION 

Following previous studies (e.g., van Dokkum et al. 2010; 
Papovich et al. 2011; Patel et al. 2013a; Leja, van Dokkum, 
& Franx 2013) we link progenitor and descendant galaxies 
by requiring that they have the same (cumulative) co-moving 
number density. Effectively, galaxies are ranked according to 
their stellar mass and we study galaxies at high redshift that 
have the same rank order as the Milky Way does at z = 0. The 
impUcit assumption is that rank order is conserved through 
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Figure 1. {a) Stellar mass density of the Universe as a function of galaxy mass, as deteiTnined from the SDSS-GALEX z = 0.1 mass function 
of Moustakas et al. (2013). (b) Evolution of the cumulative galaxy mass function from z = 0.1 to z = 3.5 (SDSS-GALEX and Marchesini et 
al. 2009). The horizontal line indicates a constant cumulative co-moving number density of 1.1 x 10"^ Mpc"'. {c) Mass evolution at a constant 
number density of 1.1 x 10" 



'^-^Mpc-^ 



cosmic time, or that processes that break the rank order do 
not have a strong effect on the average measured properties. 
As shown in Leja et al. (2013) the method recovers the true 
mass evolution of galaxies remarkably well in simulations that 
include merging, quenching, and scatter in the growth rates of 
galaxies. 

The present-day stellar mass of the Milky Way is approx- 
imately 5 X IO'^Mq (Flynn et al. 2006; McMillan 2011). 
Using the SDSS-GALEX stellar galaxy mass function of 
Moustakas et al. (2013) we find that galaxies with masses 
> 5 X 10'" M0 have a number density of 1.1 x lO'-'Mpc"-'. 
We then trace the progenitors of these galaxies by identify- 
ing, at each redshift, the mass for which the cumulative num- 
ber density is 1.1 x 10~^Mpc"^ (see Fig. [T]^). We used the 
Marchesini et al. (2009) mass functions as they are complete 
in the relevant mass and redshift range; we verified that the 
results are similar when other mass functions are used (llbert 
et al. 2013; Muzzin et al. 2013). 

The stellar mass evolution for galaxies with the rank order 
of the Milky Way is shown in Fig. [T]:. The evolution is rapid 
from z ~ 2.5 to z ^ 1 and relatively slow afterward. We there- 
fore approximate the evolution with a quadratic function of 
the form 



log(MMw)= 10.7-0.045z-0.13z 



(1) 



Based on the variation between mass functions of different 
authors, and the results of Leja et al. (2013), we estimate that 
the uncertainty in the evolution out to z ^ 2.5 is approximately 
0.2 dexFH More than half of the present-day mass was assem- 
bled in the 3 Gyr period between z = 2.5 and z = 1, and as we 
show later the mass growth is likely dominated by star forma- 
tion at all redshifts. The mass evolution is significantly faster 
than that of more massive galaxies (van Dokkum et al. 2010; 
Patel et al. 2013a), consistent with recent results of Muzzin 
etal. (2013). 

3. MILKY WAY PROGENITORS FROM Z = TO Z = 2.5 

3.1. Rest-frame Images 

Having determined the stellar mass evolution with redshift, 
we can now select galaxies in mass bins centered on this 
evolving mass and study how their properties changed. We 
selected galaxies in GOODS-North and GOODS-South as 
these fields have multi-band ACS and WFC3 imaging (from 
the GOODS and CANDELS surveys respectively; Giavalisco 
et al. 2004; Grogin et al. 2011; Koekemoer et al. 2011), as 
well as WFC3 G141 grism spectra from the 3D-HST program 

' ' We verified that changing the evolution does not affect the key results 
of this paper. 
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Figure 2. Examples of galaxies with the number density of the Milky Way at < z < 2.75. Galaxies at z « 0.015 are from the SDSS; galaxies 
at higher redshift are from the 3D-HST and CANDELS surveys. The color images were created from data in the same rest-frame bands (u and 
g) at all redshifts and have a common physical scale. Their intensities are scaled so they are proportional to mass, indicated in the top panel. 
Galaxies at high redshift have relatively low surface densities; their centers and outer parts seem to build up at the same time, at least until 

z ~ 1. 



(Brammeret al. 2012). Redshifts, stellar masses, and star for- 
mation rates were determined from deep photometric catalogs 
in these fields, combined with the grism spectra (see Bram- 
mer et al. 2012 and references therein, and R. Skelton et al., 
in preparation). The 3D-HST v2.1 catalos are w 100% com- 
plete in the relevant mass and redshift range, but we note that 
we rely almost exclusively on photometric redshifts (rather 
than grism redshifts) at z > 1.3. 

There are 361 galaxies at 0.25 < z < 2.75 in the catalogs 
whose mass is within ±0.1 dex of Mmw(z)- Images of a ran- 
dom subset of 90 are shown in Fig. |2] The images have the 
same physical scale and represent the same rest-frame filters 
(m and g). Their brightness is scaled in such a way that their 
total (m + g) flux is proportional to Mmw(z)- The rest-frame 
u and g images were created by interpolating the two ACS 
and/or WFC3 images (smoothed to the i/ieo resolution) whose 



central wavelengths are closest to the redshifted u and g filters. 

Also shown are nearby galaxies from the Sloan Digital Sky 
Survey (SDSS). We selected 40 galaxies with 0.013 < z < 
0.017 and 10.62 < logM < 10.78 from the DR7 MPA-JHU 
catalogo (Brinchmann et al. 2004), and degraded their u and 
g images to the same spatial resolution as the high redshift 
galaxies. A random subset of 10 are shown in Fig.|2] 

It is clear from Fig. |2] that present-day galaxies with the 
mass of the Milky Way have changed over cosmic time. The 
most obvious change is that galaxies became redder with time, 
particularly after z ~ 1, indicative of a decrease in the specific 
star formation rate. The galaxies also appear brighter at lower 
redshift in Fig. |2] reflecting the mass evolution of Eq. [T] A 
striking aspect of this change in brightness, and the central 
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Figure 3. Surface density profiles from z = 2.5 to z = 0, as measured from averaged, PSF-corrected rest-frame g band images in each redsiiift 
bin. The horizontal axis is linear in (a) and logarithmic in (b). The galaxy image is randomly chosen from our SDSS sample to illustrate the 
radial extent of the profiles. The main evolution is in normalization, which is determined by Mmw(z) (Eq. 1). The profile shapes are very similar 
from z ''^ 2.5 to z ~ 1, which implies that the galaxies are building up their mass at all radii. After z '--^ 1 the central regions stop growing but 
the disk continues to build. 



result of this paper, is that the bulges appear to change just as 
much as the disks, particularly at z > 1 . We do not see high 
density "naked bulges" at z ^ 2 around which disks gradually 
assembled. Instead, the central densities at z ^ 2 were much 
lower than the central densities at z ^ 0. We quantify this 
result in the remainder of the paper. 

3.2. Evolution of Surface Density Profiles 

We first analyze the surface density profiles of the galax- 
ies to study their mass growth as a function of radial distance 
from their centers. Following van Dokkum et al. (2010) we 
measured the profiles from stacked images to increase the 
signal-to-noise ratio. The galaxies were grouped in six bins 
with mean redshifts 0.015, 0.60, 1.0, 1.5, 2.0, and 2.4. Each 
bin contains at least 40 galaxies. The rest-frame u and g band 
images in each bin were normalized and stacked, aggressively 
masking all neighboring objects. 

The image stacks were corrected for the effects of the 
point spread function (PSF) following the method outlined 
in Szomoru et al. (2010). First, a two-dimensional Sersic 
(1968) model, convolved with the PSF, was fit to the stacks 
using the GALFIT code (Peng et al. 2010). Then the residu- 
als of this fit were added to the Mnconvolved Sersic model. As 
shown in Szomoru et al. (2010) this method reconstructs the 
true flux distribution with high fidelity, even for galaxies that 
are poorly fit by Sersic profiles. The resulting radial surface 
density profiles are shown in Fig. |3] The profiles are derived 
from the rest-frame g band images and scaled such that the 
total mass within a diameter of 50 kpc is equal to Mmw(z)- 
We note here that the u-g color gradients of the stacks are 
small (?s 0.1 dex"') at all redshifts, consistent with other stud- 
ies (e.g., Szomoru et al. 2013). 

There is strong evolution in the overall normaUzation of the 



profiles from z = 2.5 to z = 1 and less evolution thereafter, 
reflecting the mass evolution of Eq. [T] The evolution from 
z = 2.5 to z = 1 is strikingly uniform: the profiles are roughly 
parallel to one another in Fig. |3i), and rather than mostly 
assembling inside-out the galaxies increase their mass at all 
radii at approximately the same rate. This is in marked con- 
trast to more massive galaxies, which form their cores early 
and exclusively build up their outer parts over this redshift 
range (see Fig. 6 in van Dokkum et al. 2010 and Fig. 6 in Pa- 
tel et al. 2013a). After z ^ 1 the evolution in the central parts 
slows down but the outer parts continue to build up, consistent 
with the visual impression that around this time the classi- 
cal "quiescent bulge and star forming disk" structure of spiral 
galaxies was established (see Fig.|2]i. 

3.3. Mass Growth at Different Radii 

We explicitly show the mass growth at different radii in 
Fig. 11^. From z ~ 2.5 to z ^ 1 the mass within r = 2 kpc 
and the mass outside of 2 kpc increase at approximately the 
same rate, mirroring the increase in the total mass over this 
time interval. Again, this trend is qualitatively different from 
that seen in more massive galaxies: after z '--^ 2 the mass 
within 2 kpc is constant to within 0. 1 dex for galaxies with 
log(M/MQ)(z = 0) = 11.2 (see Fig. 7 of Patel et al. 2013a). 
After z ~ 1 the trend is similar to that for more massive galax- 
ies: the central mass is approximately constant, whereas the 
mass at large radii continues to grow. 

In Fig.|4j3 we express the growth in mass as an (implied) star 
formation rate. The star formation rate was calculated directly 
from Eq. [T] with a x 1.35 upwards correction to account for 
mass loss in windsQ The implied star formation rate is ap- 

" This factor is the mass loss after 2 Gyr for a Kroupa (2001) IMF. 
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Figure 4. (a) Comparison of the mass growth in the central regions to the growth at larger radii. The galaxies grow at all radii until z ~ 1, after 
which the mass inside r = 2 kpc remains roughly constant, (b) Implied evolution of the star formation rate. Data points are the mean measured 
star formation rates of the galaxies in each redshlft bin, from the 3D-HST v2. 1 catalogs (Skelton et al., in preparation). There is an excellent 
match between the black curve and the points, indicating that mergers are not required to explain the mass evolution of large spiral galaxies. 



proximately constant at 10- ISM© yr 

and then decreases rapidly to < 2 Mq yr 

of this star formation history is well approximated by 



' from z ~ 2.5 to z '--^ 1 
' at z = 0. The form 



log(l + SFR) = 0.26 + 0.92z-0.23z^ 



(2) 



This evolution is qualitatively consistent with the color evolu- 
tion seen in Fig. |2l as the galaxies become rapidly redder in 
rest-frame m-^ at z < 1. More to the point, Eq.|2]is consistent 
with the actual star formation rates of the galaxies: the points 
with errorbars in Fig. |4]3 show the mean star formation rates 
of the galaxies that went into the analysis, as obtained from 
SED fits (see Kriek et al. 2009, and Skelton et al., in prepara- 
tion). This consistency implies that the assembly history can 
be fully explained by star formation and that mergers likely 
played a minor role. This can, again, be contrasted with more 
massive galaxies, as star formation is not sufficient to explain 
their growth after z ^ 1 .5 (van Dokkum et al. 2010). 

3.4. Structural Evolution 

Finally, we quantify the implications of our results for the 
structural evolution of galaxies with the present-day mass of 
the Milky Way. As the mass growth is mostly independent 
of radius, we expect the structure of the galaxies to remain 
more or less the same over cosmic time. The evolution of 
the GALFIT-derived structural parameters of the stacks (see 
§ 3.2) is shown in Fig.|5] 

The effective radii and Sersic indices have indeed changed 
relatively little since z ^ 2.5, particularly when it is consid- 
ered that the galaxies increased in mass by a factor of '^ 10 
over this time. The radius increased by a factor ^1.8 and 
the Sersic index changed from « ~ 1.5 to « ~ 2.5. The red 
curves show the change in these same parameters for high 
mass galaxies, calculated in the same way (Patel et al. 2013a). 
Even though the progenitors of today's massive galaxies in- 
creased their mass by only a factor of ~ 3 over this redshift 
range they show much more dramatic structural evolution. 

This point is emphasized in Fig. |5};,d which compares the 
structural evolution to the mass evolution for both classes of 



galaxies. The sizes of massive galaxies grow as r^ ex M (van 
Dokkum et al. 2010; Hilz et al. 2013; Patel et al. 2013a), 
whereas those of galaxies with the mass of the Milky Way 
grow as re oc M°^. This slope is similar to that of the size- 
mass relation of late type galaxies (e.g., Shen et al. 2003). We 
note that an increase in Sersic index does not imply growth 
of a classical bulge for either class of galaxy (see also Nelson 
etal. 2013). 



4. DISCUSSION 

In this Letter we have demonstrated that it is possible to ob- 
tain a description of the formation of galaxies with the mass 
of the Milky Way all the way from z '^ 2.5 to the present. We 
find that these galaxies built up ^ 90 % of their stellar mass 
since z ^ 2.5. The build-up can be fully explained by the mea- 
sured star formation rates of the galaxies, and does not require 
significant merging. The key result of our paper is that most of 
the mass growth took place in a uniform way, with the galax- 
ies increasing their mass at approximately the same rate at 
all radii. Our results are therefore inconsistent with models 
in which bulges are fully assembled before disks: we do not 
find "naked bulges" at high redshift. Instead, they are con- 
sistent with models in which bulges (and presumably black 
holes) were largely built up at the same time as disks, through 
short-lived peaks in the accretion rate, bar instabilities, clump 
migration, or other processes (e.g., Debattista et al. 2006; 
Dekel et al. 2009). The implied star formation rate declines 
precipitously after z ^ 1 , particularly in the central w 2 kpc of 
the galaxies. By z = we are left with quiescent bulges and 
slowly star forming disks. 

Many other studies have reached similar conclusions us- 
ing independent arguments; here we limit the discussion to 
a handful examples. Wuyts et al. (2011) and Nelson et al. 
(2013) find that star formation at high redshift typically oc- 
curs in disks. Nelson et al. (2012) find that galaxies begin 
to build inside-out at z ^ 1. As noted in § 1 Genzel et al. 
(2008), Forster Schreiber et al. (201 1), and others have iden- 
tified thick, clumpy star forming disks at z ^ 2. Finally, the in- 
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Figure 5. Effective radius and Sersic index as a function of redshift and mass, for Milky Way progenitors (blue) and more massive galaxies (red, 
taken from Patel et al. 2013a). Galaxies like the Milky Way have undergone much less structural evolution than the giant elliptical galaxies that 
populate the high mass end of the mass function. 



ferred star formation history (Eq. 2) is consistent with the one 
derived by Behroozi, Wechsler, & Conroy (2012) for 10 '^Mq 
halos, but is more peaked than the history inferred by Moster, 
Naab,& White (2013). 

It is tempting to compare our results directly to known prop- 
erties of the Milky Way itself; e.g., Eq.|2]implies a z = star 
formation rate of ^ 1 Mq yr"', in reasonably good agreement 
with that of the Milky Way (Robitaille & Whitney 2010). We 
note, however, that the Milky Way has a relatively low bulge- 
to-disk ratio for its mass (e.g., McMillan 201 1). Furthermore, 
the Milky Way, like any other galaxy, has had a unique his- 
tory and it is fundamentally hazardous to apply the statistical 
analysis of samples of distant galaxy samples to an individual 
nearby galaxy (see, e.g.. Fig. 1 of Leja et al. 2013). 

As noted in previous Sections, the formation process of 
galaxies with logM ?» 10.7 appears to be very different from 
that of more massive galaxies. Massive galaxies formed 
inside-out, with their extended wings assembling after forma- 
tion of a compact core at high redshift. It will be interesting 
to see if galaxy formation models can reproduce both types 
of behavior seen in Fig.|5j e.g., it may be that (minor) merg- 



ers lead to inside-out growth whereas gas accretion leads to 
uniform growth. 

This study can be extended and improved in many ways. 
Most importantly, we have largely ignored systematic uncer- 
tainties in our analysis. Among the uncertainties are the low 
mass end of the mass function at z > 2 (see, e.g., Reddy & 
Steidel 2009); possible errors in the number density selection 
technique (Leja et al. 2013); systematic errors in redshifts 
and/or masses in the 3D-HST v2.1 catalogs; and the conver- 
sion of light-weighted to mass-weighted profiles. We have 
also ignored the spread in galaxy properties at fixed mass (see, 
e.g., Baldry et al. 2006; Franx et al. 2008, and Fig. 2). Fi- 
nally, our analysis is, by its nature, indirect: we do not actually 
observe the formation of different parts of the galaxies but in- 
fer this from changes in their stellar surface densities. Stellar 
migration and other processes almost certainly altered the or- 
bits of stars after their formation (Roskar et al. 2008). Deep, 
direct observations of spatially-resolved gas distributions at 
high redshift, particulary in the crucial epoch 1 < z < 2.5, are 
needed to disentangle formation and migration, and to shed 
light on the physical processes that are at work (e.g.. Nelson 
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et al. 2012, 2013; Freundlich et al. 2013). 
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